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ABSTRACT 

Confexf. 

Aims. We develop a model of time delays between the continuum bands in the Narrow Line Seyfert 1 galaxy Mrk 335 to explain the observed 
delays measured in this source. 

Methods. We consider two geometries: an accretion disk with fully ionized warm absorber of considerable optical depth, located close to the 
symmetry axis, and an accretion disk with a hot corona. Both media lead to significant disk irradiation but the disk/corona geometry gives 
lower values of the time delays. 

Results. Only the disk/corona models give results consistent with measurements of Sergeev et al., and a low value of the disk inclination is 
favored. The presence of an optically thick, fully ionized outflow is ruled out at the 2-cr level. 
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1. Introduction 

Active galactic nuclei (AGN) are powered by accretion onto 
a central massive black hole, and the accretion is believed to 
proceed predominantly through the accretion disk. However, 
the shapes of the broad band spectra of most AGN are much 
broader than predicted by the classical disk model of Shakura 
& Sunyaev (1973). The deviation is not as strong in the case 
of bright quasars (e.g. Zheng et al. 1997, Laor et al. 1997, 
Koratkar & Blaes 1999, Czerny et al. 2004, Shang et al. 2005), 
but in the case of Seyfert 1 galaxies the role of the X-ray emis- 
sion is essential and the presence of the accretion disk is far less 
obvious when looking at the optical data. The spectral slopes 
are much redder than the standard slope of p = 1/3 of the 
Shakura & Sunyaev (1973) disk in the F v oc v p representation 
(e.g. Edelson & Malkan 1986; Kong et al. 2004). The analy- 
sis of the X-ray spectral features may support the view that the 
disk is disrupted in the inner part of the flow (a relatively nar- 
row iron Ka is line preferred, e.g. Yaqoob et al. 2001, Reeves 
et al. 2004, although a weak broad component cannot be ruled 
out with high significance, e.g. Yaqoob & Padmanabhan 2004, 
Yaqoob et al. 2005). 

Narrow Line Seyfert 1 galaxies (NLS1) are an intermediate 
type of objects: their spectra are dominated by disk emission, 
as in bright quasars (see e.g. PG121 1 + 143, Czerny & Elvis 
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1987; RE J1034+396, Puchnarewicz et al. 2001; Mrk 335, 
Pounds et al. 1987), and the disk is likely to extend down to the 
marginally stable orbit (as inferred from the iron line shape in 
MCG -6-30-15, Tanaka et al. 1995, Fabian et al. 2002; 1H0707- 
495, Fabian et al. 2004; PG121 1 + 143, Janiuk et al. 2001, IRAS 
13349+2438, Longinotti et al. 2003) but the deviations from 
the standard disk are quite strong as well. NLS 1 galaxies usu- 
ally have profound soft X-ray excesses and/or steep (soft) X- 
ray spectra (e.g. Brandt et al. 1997, Leighly 1999, Romano et 
al. 2004, Boiler 2004), and their X-ray spectra frequently show 
warm absorber features (e.g. in MCG-6-30-15, Turner et al. 
2004; NGC 4051, Collinge et al. 2001, Ton S 180, Rozariska 
et al. 2004; PG 1211 + 143, Pounds et al. 2003). It is therefore 
interesting to explore which additional component, apart from 
an accretion disk, might be responsible for the overall spectral 
shape. 

We select one of the most studied NLS1 galaxies, Mrk 
335 (z = 0.026), for the modeling. This is a typical NLS1 
object, with a narrow H/3 line (1640 km s , Wang et al. 
1996), and variable Big Blue Bump with a large soft X-ray 
excess (Edelson & Malkan 1986, Pounds et al. 1987, Turner 
& Pounds 1988, Turner et al. 1993, Nandra & Pounds 1994, 
Leighly 1999), possibly a relativistically smeared iron Ka line 
(Gondoin et al. 2002, Crummy et al. 2006), and direct evidence 
of a weak partially ionized warm absorber (Turner et al. 1993, 
Reynolds 1997, Bianchi et al. 2001). Since the optical/UV 
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slope is roughly a power law (F v cc v° , with a = -0.64 + 0.01, 
Constantin & Shields 2003) and does not constrain the model 
parameters very strongly, we concentrate here on modeling the 
observed time delays between the UBVRI bands measured for 
this object. The source was monitored for several years at the 
Crimean observatory and the time delays between the UBVRI 
bands were measured (Doroshenko et al. 2005, Sergeev et al. 
2005). We describe the model in Sect. [2] In Sect. [3] we describe 
its application to Mrk 335, and in Sect. [4] we discuss the con- 
sistency of the results with independent estimates of the central 
black hole mass and the optical depth of the warm absorber. 

2. Model 

The UV/optical/IR continuum of Mrk 335 varies significantly 
(Sun et al. 1994, Kassebaum et al. 1997, Peterson et al. 1998, 
Tao et al. 2004, Klimek et al. 2004). The measured delays of the 
continuum between the U band and BVRI bands increase with 
the wavelength difference (Sergeev et al. 2005, Doroshenko 
et al. 2005). Such variations in AGN spectra are most likely 
caused by the irradiation of the outer parts of the accretion 
disk responsible for the optical/UV emission by the EUV and 
X-ray radiation generated in the innermost part of the disk, 
as discussed by Collin-Souffrin (1991) (see also Courvoisier 
& Clavel 1991, Krolik et al. 1991, Collier et al. 1998 and 
Oknyanskij et al. 2003). The direct irradiation of the outer parts 
of the disk by the inner disk is unlikely (Loska et al. 2004). The 
reason is that AGN disks are geometrically very thin, much 
thinner than GBH disks. In Seyfert 1 galaxies the irradiation is 
possible due to the extended optically thin X-ray emitting in- 
ner region, which increases the geometrical factor (Rokaki et 
al. 1993, Chiang 2002). However, in NLS1 the standard, op- 
tically thick disk extends to the marginally stable orbit, rather 
than being disrupted in its inner part, and therefore the presence 
of an optically thin inner flow is rather unlikely. 

Therefore, we suggest here that additional material may ef- 
ficiently redirect the nuclear emission towards the outer disk. 
We consider in detail two geometries: a warm absorber type 
medium, predominantly occupying the region around the sym- 
metry axis, and the disk wind/corona, located closer to the disk 
plane. 

2.1. Disk with a warm absorber 

We apply here the simple model of disk irradiation described 
in Czerny et al. (2003). We neglect the true absorption effect 
due to the warm absorber since (i) we are mostly interested in 
the fully ionized fraction of this material (ii) the absorption is 
limited to the 0.5 - 2 keV energy band, and the scattering ef- 
fect (working in the whole energy band) is more effective than 
absorption and reemission in redirecting outgoing photons to- 
wards the disk. We assume that most of the scattering medium 
is located near the symmetry axis in the form of a cone (see 
Fig. [TJ. The warm absorber parameters are: the fraction of the 
sky covered by the cone, f mt (equivalently, the opening angle of 
the cone), the smallest and the largest distance, z m m and Zmax, 
of the scattering material measured along the symmetry axis, 
the parameter 6 that measures the density distribution of the 






Fig. 1. The two geometries used: a disk/warm absorber geom- 
etry (panel a) and disk/corona geometry (panel b). Differences 
in the location of the scattering medium lead to differences in 
the time delays between the continuum bands. 

warm absorber ( p(z) °c z~ s ), and the total optical depth of the 
medium for scattering, r. 

The local irradiation flux is given by 



F(r) 



£ 



T(r,z)dz, 



where 
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with L being the total luminosity of the source. This approach 
is much less time consuming than the more accurate modeling 
done by Loska et al. (2004). However, our model is still capa- 
ble of producing reliable results. In the most extreme case of 
irradiation considered for the source RE J 1034+396, the sim- 
plified numerical scheme underestimates the flux in the optical 
range by 25 % (see Fig. [2]) but the overall shape is preserved; in 
moderate illumination cases, with higher 5, the error does not 
exceed a few percent. 

The disk internal emission is described in the Newtonian 
approximation, and parameterized by the black hole mass, M, 
and the dimensionless accretion rate, m, defined with the ef- 
ficiency factor of 1/12 included (appropriate for Newtonian 
boundary condition, Shakura & Sunyaev 1973), i.e. 
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The illuminated disk is assumed to radiate locally as a black 
body. This simplifying assumption allowed us to save computer 
time, and is fully justified. However more detailed models of 
the disk emission do show a departure of the local spectrum 
from a black body, but this is only if the opacities of heavy 
elements are neglected. As shown e.g. in Hubeny et al. (2001) 
and Madej & Rozariska (2000), taking these opacity effects into 
account minimizes the difference with respect to a blackbody 
spectrum (apart from the soft X-ray part). 
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Fig. 2. Comparison of the full 3-D irradiation computations 
with self-irradiation included (Loska et al. 2004, thin contin- 
uous line) with the results of the simplified scheme used here 
(thick continuous line) in the case of extreme illumination re- 
quired by the source RE J 1034. The bare disk model is marked 
with a dotted line. Parameters: M - 6.3 x 10 5 M o , m = 0.84, 

5 = 0.0, T = 0.6, Zmin = lORschw, Zmax = 1000R S clm; fwa = 0.5. 



The emitted disk spectrum passes through the warm ab- 
sorber, so the observed spectrum is reduced by the factor 
exp(-r) with respect to the intrinsic spectrum. In this way we 
neglect the possible dumpiness of the warm absorber, which 
in principle may decouple the specific line of sight from the 
average scattering probability. 

The time delays are calculated assuming that, for a given 
radius, all irradiating flux comes from a single representative 
distance on the symmetry axis. This distance is calculated as 
the effective distance, z e ff in the following way 



Zeff 



\ Zm ""T(r,z)zdz 

''Zmin 

f z Zmax T(r,z)dz' 



(4) 



The delay of the emission from a given radius with respect to 
the emission from the innermost part of the disk is calculated 

as 



2U/2 



Af(r) = r sin i + z eff + {zf eff + r z ) 



2.2. Disk with a flaring corona 



(5) 



The accretion disk is also likely to develop an extended corona 
or a wind-type outflow. In this case the additional hot mate- 
rial is located predominantly close to the disk surface but high 
enough to intercept some of the nuclear flux much more effec- 
tively than the disk surface itself. There may be some dissipa- 
tion taking place directly in the wind-corona. Both effects will 



increase the local disk temperature. Since there are no well jus- 
tified models of this situation, we model such a disk assuming 
that a fraction of the disk, between r mi „ and r max , is irradiated 
by additional flux. This irradiating flux is simply parameterized 
as a power law 



(6) 



where the normalization coefficient A is introduced through the 
parameter / describing the ratio of the total energy in this com- 
ponent to the total energy due to accretion 
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Again, we tested the analytical approach against the 3-D nu- 
merical computations of the scattering in a corona which has 
a density profile in the radial direction but is uniform in the 8 
direction (in spherical coordinates) and has a constant open- 
ing angle. The results are illustrated in Fig. [3] Our analytical 
model is approximately correct although the exact results pre- 
dict some irradiation effect in the innermost part of the disk 
(below r m i n , roughly constant with radius) as well as above r max 
(decreasing roughly as r~ 3 ). The irradiation is significantly dif- 
ferent from a point-like source since in our model we have an 
extended region with incident flux decreasing as an arbitrary 
power of radius (physically determined by the density distribu- 
tion in the corona) while in the case of a point-like irradiation 
the incident flux is always constant below a certain radius and 
decreases as r~ 3 above it. 

In computations of time delays we neglect the extension of 
the coronal medium in the vertical direction, so the emission 
from a given radius is delayed with respect to the central re- 
gions by 



Af(r) = r(sin/ +1). 



2.3. Time delays in UBVRI bands 



(8) 



In order to determine the time delay at a given wavelength with 
respect to the emission from the innermost part of the disk 
we have to connect a given radius with a representative wave- 
length. In the simplest approximation, this can be done taking 
into account that a local black body emission with the temper- 
ature T peaks at 2.1kT in a log F Y vs. log v diagram. However, 
if the spectrum is not of a single black body type, the domi- 
nant contribution may not be well determined by such a crite- 
rion (for example, for a power law spectrum a factor before kT 
depends on the spectral slope; see Siemiginowska & Czerny 
1989). Since the shapes of the model spectra are more complex 
than a single power law, we use a numerical approach. In order 
to connect a given radius with a representative wavelength, we 
calculate a model with a narrow gap in the disk around a cer- 
tain radius, then we compare the spectrum to the model with- 
out a gap. We find the wavelength where the ratio of the gap 
model spectrum to the full model spectrum was the lowest. By 
repeating such a procedure we identify a set of radii with a cor- 
responding sequence of wavelengths. 
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Fig. 3. Comparison of the 3-D irradiation computations of scat- 
tering by the corona with radial density p oc r~ 08 , extending 
from - l30Rschw to r out = 10 Rschw, total optical depth in 
the radial direction 0.734 and constant opening angle of 15° 
measured from the equatorial plane (dashed line) with the re- 
sults of the simplified scheme used in the present paper (con- 
tinuous line) for equivalent parameters, i.e. the same r in , r out , 
black hole mass and accretion rate, M = 10 7 M Q , m = 0.95, 
but with y — 1.8 and / = 0.095. Such a value of the factor / 
is equivalent to the scattering probability given by the optical 
depth and the azimuthal extension of the corona in numerical 
computations. It can be shown analytically that the power law 
distribution of the coronal density with a given slope leads to 
a disk region irradiated by the flux decreasing as a power law 
with index larger by 1 . 

2.4. Data 

The Galactic hydrogen column in the direction of Mrk 335 is 
N H = 3.8 x 10 20 cnr 2 (Veron-Cetty et al. 2001). The reddening 
can therefore be estimated by taking the standard relation for 
the interstellar medium in the Galaxy N H /E(B-V) - 5.8 xlO 21 
(Bohlin et al. 1978) which gives E(B - V) = 0.066 and the 
data were corrected accordingly. The continuum fluxes used 
for model fitting are given in Table Q] 

The source was systematically monitored in Crimea for 10 
years (Doroshenko et al. 2005, Sergeev et al. 2005). The mea- 
sured delays between the various optical bands are given in 
Table El 

The measured delays generally depend on the method (peak 
or centroid), and the errors in measurements with respect to the 
U band are much larger than those measured with respect to 
the B band. The I vs. B delay of Sergeev et al. (2005) gives an 
upper 1 cr limit of 2.5 - 2.6 day, independent of the method. 

The delays measured by Doroshenko et al. (2005) are sys- 
tematically larger. All models are consistent with 1 <x upper 
limits of more than 17 days of delay between the U and I band. 



Table 1. Continuum data points of Mrk 335 used for model 
fitting 



logv 


flux [mjy] 


U 


14.23 


B 


13.99 


V 


15.51 


R 


18.09 


I 


19.49 


Jl 


21.00 


J2 


22.10 



The data points are taken from NED (McAlary et al. 1983, Spinoglio 
et al. 1995). 
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Fig. 4. Representation of the optical data by two bare disk 
models: M = 7 x 10 8 M o , m = 0.0015 (dotted line), and 
M — 5 x 10 7 M Q , m — 0.1 (long dash-dot line), an example of 
the fit with a warm absorber: M = 1.78 X 10 7 M o , m = 0.6, 

8 = 1.4, T = 0.8, Zmi„ = 50R Sc hw, Zmax = 40000R Schw , 

cos i = 1 (continuous line), and an example of the fit with 
corona: M = 1.0 x 1O 7 M , m = 0.95, y = 1.8, / = 0.095, 
r min = l30R Schw , r max = 10000# Sc fcw. cos/ = 1 (long dashed 
line). The solid squares mark the data points (cf. Table[T). 

However, we consider this delay measurements as much less 
certain than the delays given by Sergeev et al. (2005) since the 
measurement of U flux is much less certain than the measure- 
ment of B, when a relatively small telescope located in a poor 
astronomical site is used. 

3. Results 

3.1. Bare disk model 

The presence of time delays in Mrk 335 between the optical 
bands, increasing with the wavelength separation, clearly sug- 
gests the role of reprocessing in the formation of the optical 



B. Czerny and A. Janiuk: Modeling the UBVRI time delays in Mrk 335 
Table 2. Measured time delays in Mrk 335 



5 



Band Doroshenko et al. 


Doroshenko et al. 


Sergeev et al. 


Sergeev et al. 


peak 


centroid 


peak 


centroid 


U 








B 0.4!" 

v e.o:ii 

Rc 5.0+lf 
I/Ic 7.6 + !', 2 


9 9+6.8 

4.5*« 
8.3^' 
7.6+%> 


43 +0 - 83 

-0.81 
■ 6 * JU -1.76 

2.33!?;» 


0.86^* 

1 4^+0.94 
1- ^ J -1.26 
1 74+0.78 
1 - '^-1.25 



The delays measured with respect to the U band in Doroshenko et al. (2005) and with respect to the B band in Sergeev et al. (2005). Either the 
standard or Cousins filter system was used (Rc,Ic). Quoted measurement uncertainties are 1 a errors. 



spectra. However, as a reference we consider a simple fit to the 
spectrum by a non-irradiated standard disk. 

The spectrum is relatively blue at longer wavelengths (p * 
-0.2), but it is still not as steep as expected from the standard 
disk, and it flattens at shorter wavelengths tops -0.5. Such 
a curved spectrum can be fitted by a standard disk if the black 
hole mass is high, so the maximum disk temperature is within 
the optical/UV band. An example of the data fit without the 
presence of the warm absorber is shown in Fig. [4] (dotted line). 
The model does not fit the near-IR data points but this can be 
explained by starlight contamination. However, such a model 
cannot account for the total luminosity of the source (the model 
bolometric luminosity is only 1.3xl0 44 erg s _1 while the source 
bolometric luminosity was estimated to be much higher, 1.3 X 
10 45 erg s- 1 ; Pounds et al. 1987, or 6 x 10 45 erg s _1 ; Edelson & 
Malkan 1986). The required mass is much higher than the value 
obtained from the reverberation (M = 1.42 + 0.37 x 10 7 M o , 
Peterson et al. 2004). 

On the other hand, if we fix the accretion rate at some larger 
value, as expected for NLS 1 galaxies, and we adjust the black 
hole mass to the required emission level, the resulting black 
hole mass is not so large (e.g. M = 5 x 1O 7 M for adopted m = 
0.1) but the predicted spectrum is much too blue to represent 
the data points (see the steepest line in Fig.|4|i. 

Therefore, the bare disk model is unlikely to represent well 
the Mrk 335 data. Most of the older models based on the fit to 
optical/UV data alone would now face the same problem, as the 
masses inferred by such models are high (e.g. Sun & Malkan 
1989,M = 8.9xlO 7 M fora top view, and it is largerfor higher 
inclination angles; Zheng et al. 1995, M = 5 x 1O 7 M ). 

3.2. Irradiated disk model 

Both geometries, irradiation due to the scattering by a fully ion- 
ized warm absorber and due to scattering by the disk corona, 
can account for the overall shape of the spectrum. We did not 
try to differentiate between the models on the basis of the qual- 
ity of the fit because the observational points can have some 
systematic errors due to slight internal reddening, starlight con- 
tamination and the contribution of the small blue bump (i.e. Fe 
II and Balmer continuum) to the continuum. Also the theoreti- 
cal predictions are not precise; various treatments of the radia- 
tive transfer in the disk surface can lead to differences between 
the models of the order of 20% in the U band if the normal- 
ization in the I band is fixed (see e.g. Ross et al. 1992, Hubeny 



et al. 2001). Since for each mass, accretion rate and inclination 
angle we fit other model parameters, we obtain a satisfactory 
data representation for a broad range of these basic parame- 
ters. The results for the time delays for acceptable values of 
mass and accretion rates are given in Table [3] and Table H] for 
the two models: the warm absorber and disk/corona system, 
respectively. 

In both cases the acceptable solutions create an extended 
elongated region in M - m parameter space since an increase 
in the black hole mass can be compensated by the decrease in 
the accretion rate in order to normalise the spectrum consis- 
tently with the data (see for example the discussion by Tripp 
et al. 1994 of the bare disk data fitting). The location of this 
region depends on the inclination angle: larger disk inclina- 
tions allow for larger masses for a given accretion rate. The 
width of this elongated region is provided by the additional 
model parameters, and it is much larger in the disk/warm ab- 
sorber scenario than in the disk/corona scenario. For example, 
if we fix the black mass at 1O 7 M and the inclination angle at 
cos i = 1 .0, the acceptable fits in the disk/warm absorber ge- 
ometry are obtained for m = 0.6 (5 — 1.4, r = 0.25, f wa = 1.0, 
Zmin = l30R Schw , Zmax = 40000/? Sc/w ,), for m = 0.8 (6 = 1.4, 
r = 0.30, f wa = 0.6, Zmin = UQRschw, Zmax = l0000R Sc hw), and 
for m = 0.8 (6 = 1.6, r = 0.15, f wa = 0.6, Zmm = 310R Schw , 
Zmax = 4OQ0QRschw)- In the disk/corona geometry, the accept- 
able fits are for m = 0.8 (y = 1.8, / = 0.13, r min - 5QRs c hw, 
r max = 10000/? Sc/n ,) and for m = 1.0 (y = 1.8, / = 0.095, 
r m in = l30R Sc hw, r ma x = l0000R Sc hw)- 

We give in Table [3] and Table H] the theoretical delays be- 
tween the B and I band. The delays are systematically larger 
for the disk/warm absorber geometry than for the disk/corona 
geometry since the light travel time through the warm absorber 
and back to the disk contributes to the delays. 

We compare these delays to the observations summarized 
in Table|2] Only predictions of the disk/corona geometry at low 
inclination (cos / = 1) are consistent with the upper 1 cr limit 
of Sergeev et al (2005) These models cover a broad range of 
accretion rates (from 0.2 to 1 in dimensionless units) and black 
hole masses range from 1O 7 M to 3 x 1O 7 M . 

An example of the wavelength-dependent time delays for 
one of the disk/corona models is shown in Fig. [5] We see a 
good agreement between the model and the data. The delays 
increase with the wavelength separation in a similar way in the 
model and in the data. 
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Table 3. Time delays between B and I band in days from the 
disk/warm absorber system 



log M: 


6.5 


6.75 


7.0 


7.25 


7.5 


7.75 


8.0 








cos i 


= 1.0 








log m 


0.2 


X 


X 


X 


X 


3.0 


4.3 


X 


0.4 


X 


X 


X 


3.5 


3.6 


7.7 


X 


0.6 


X 


X 


3.7 


6.1 


5.6 


X 


X 


0.8 


X 


X 


3.2 


3.3 


6.2 


8.3 


X 


1.0 


X 


2.9 


3.0 


4.5 


X 


X 


X 








cos i 


= 0.75 








log m 


0.2 


X 


X 


X 


X 


6.0 


4.4 


9.3 


0.4 


X 


X 


X 


5.2 


9.7 


8.0 


11.0 


0.6 


X 


X 


5.1 


4.6 


8.7 


X 


X 


0.8 


X 


X 


6.1 


9.3 


7.8 


X 


X 


1.0 


X 


X 


7.0 


6.3 


9.5 


X 


X 








cos i 


= 0.50 








log m 


0.2 


X 


X 


X 


X 


X 


8.2 


12.0 


0.4 


X 


X 


X 


X 


7.3 


9.0 


13.0 


0.6 


X 


X 


X 


6.6 


10.0 


12.0 


14.0 


0.8 


X 


X 


7.6 


6.9 


13.0 


14.0 


17.0 


1.0 


X 


X 


5.9 


6.7 


17.0 


X 


X 



Accretion rate is given in dimensionless units, assuming standard ef- 
ficiency of accretion onto a non-rotating black hole. Other model pa- 
rameters were found as the best fit to the optical/IR data of Mrk 335. 
Signs 'x' mark the values of black hole masses and accretion rates for 
which no acceptable model can be fitted to the data. 



Table 4. Time delays between B and I band in days from the 
disk/corona system 



log M: 


6.5 


6.75 


7.0 


7.25 


7.5 


7.75 


8.0 








cos i = 


: 1.0 








log m 


0.2 


X 


X 


X 


X 


1.7 


X 


X 


0.4 


X 


X 


X 


2.1 


X 


X 


X 


0.6 


X 


X 


X 


1.9 


X 


X 


X 


0.8 


X 


X 


1.9 


X 


X 


X 


X 


1.0 


X 


X 


2.0 


X 


X 


X 


X 
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X 


X 


X 


X 
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X 


X 


X 


X 


3.3 


X 


X 
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X 


X 


X 
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X 


X 


X 


3.7 


X 


X 


X 


1.0 


X 


X 
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X 


X 


X 








cos i = 


0.50 








log m 


0.2 


X 


X 


X 


X 


X 


4.0 


X 


0.4 


X 


X 


X 


X 


X 


X 


X 


0.6 


X 


X 


X 


X 


5.3 


X 


X 


0.8 


X 


X 


X 


5.8 


X 


X 


X 


1.0 


X 


X 


4.5 


5.3 


X 


X 


X 



Accretion rate is given in dimensionless units, assuming standard ef- 
ficiency of accretion onto a non-rotating black hole. Other model pa- 
rameters were found as the best fit to the optical/IR data of Mkr 335. 
Signs 'x' mark the values of black hole masses and accretion rates for 
which no acceptable model can be fitted to the data. 



None of the disk/warm absorber models is acceptable 
within a 1 cr error. 

If a larger 3 cr error is allowed, the upper limit for the I vs. 
B delay increases to 7.1 days (centroid measurement) and to 
7.0 day (peak measurement). All disk/corona models, as well 
as several low inclination examples of the disk/warm absorber 
geometry, satisfy this requirement. 

Although the numerical method allows us to determine the 
time delays precisely from the point of view of the adopted 
algorithm, the actual errors of such determinations are signifi- 
cant. Using the method outlined above, we tested the radial ex- 
tension of the region responsible for emission at a given wave- 
length. For a specific model, the radial extension of the zone 
responsible for 50% of the flux in B band was from 25QRs c hw 
to 1260Rschw, and the zone responsible for 50% of the flux in 
I band lay between 560Rsdnv and 440QRs ( -/ m .. It is therefore not 
surprising that the errors of the measured time delays are also 
large. 

We also checked whether the dependence of the variabil- 
ity level on the wavelength is well reproduced by any of the 
models. We considered two luminosity states in a single rep- 
resentative model: one with the accretion rate (and luminosity) 
enhanced in the innermost 10 Rschw and one with a suppressed 
innermost accretion rate. The enhancement/suppression factor 
was adjusted to the ratio of the minimum to maximum flux 
in B band given by Doroshenko et al. (2005). We used their 



photometric measurements for that purpose. The basic model 
parameters were as in Fig.0] 

We give the results in Table [5] for the disk/corona model, 
but the results for the disk/warm absorber geometry are almost 
the same. In the data, the variability amplitude decreases to- 
wards IR but in the model it actually increases. This is due 
to the fact that we assumed a constant intrinsic disk emission 
at larger radii and variable irradiation, and that the role of ir- 
radiation increases towards the IR in order to reproduce the 
observed spectra that are much redder than the intrinsic disk 
spectrum. We can assume that a certain amount of starlight con- 
tributes to the spectrum in the VRI band. However, the amount 
of starlight required to suppress the variability is very large (see 
column 4 in Table |5). If the spectrum is corrected for such an 
amount of starlight the intrinsic spectrum becomes bluer than 
in the standard disk, leaving no space for the effect of irradia- 
tion and time delays. We cannot perform self-consistent com- 
putations for the full model, with starlight iterations (i.e adopt- 
ing a certain starlight distribution, finding new best solution 
for the spectrum, checking the variability properties, correct- 
ing the starlight again etc), since this is too time consuming. 
However, we estimated the probability of finding a consistent 
solution satisfying both the requirement of a red spectral slope 
and variability amplitude decreasing strongly with increasing 
wavelength. We used for that purpose an analytical power law 
description for all three spectral component independently, as- 
suming an arbitrary power law index for the starlight contribu- 
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Fig. 5. Time delays measured with respect to the B band. Data 
points: centroid delays from Sergeev et al. (2005). Model: 
disk/corona geometry with parameters as in Fig. [4] 

Table 5. The variability properties in the BVRI band 





F 'max I min 




starlight 




measured 


disk/corona model 


required 


B 


1.550 


1.550 





V 


1.384 


1.667 


36 


R 


1.358 


1.703 


42 


I 


1.343 


1.819 


55 



The measurements are from Doroshenko et al. (2005), photometric 
data, the parameters of the disk/corona model are as in Fig. [4j the 
fourth column gives the amount of starlight required to dilute the vari- 
ability to match the observations. 



tion, an index of 1/3 for the disk, an index of -1 for the spectral 
component due to irradiation, and we searched the parameter 
space of the starlight index and the two normalization factors, 
determining starlight to disk and irradiation to disk ratios. In 
all solutions with variability larger in B band than in I band, 
the difference in the variability level was by 1 - 2 %, the IR 
slope in the final spectrum was unacceptably steep (with power 
law index above 2.7), and the overall variability level in B band 
was too low (F max /F m i„ ~ 1.1) for models with the IR index 
smaller than 3.0. 

The variability at longer wavelengths may be smeared, and 
the amplitude effectively lowered due to the larger reprocessing 
area. More advanced time-resolved modeling is necessary to 
verify this possibility. 

4. Discussion 

We considered two geometries that can provide the required ef- 
ficient irradiation of the outer parts of an accretion disk respon- 



DISK/CORONA GEOMETRY 



top view observer 




DISK/WARM ABSORBER GEOMETRY 



Fig. 6. Schematic illustration why the time delays are system- 
atically larger in the disk/warm absorber geometry of the scat- 
tering medium than in the disk/corona geometry. 

sible for the optical emission: a disk/warm absorber geometry 
and a disk/corona geometry. In the first case the scattering ma- 
terial is located relatively far from the disk, around the symme- 
try axis, while in the second geometry the scattering material is 
located closer to the disk surface. 

We compared the predictions of these two geometries to 
the measurements of the time delays between the various con- 
tinuum bands. The expected time delays are generally shorte 
by a factor of 2-3 in the case of the disk/corona geometry than 
in the disk/warm absorber geometry. This is due to the fact that 
for a disk radius of interest, r, the effectively scattering coronal 
region is located roughly at the disk radius distance from the 
center while the most effectively scattering warm absorber re- 
gion is located at the distance r along the symmetry axis and the 
overall distance traveled by photons is roughly r + yfr — 2Ar 
before the radiation reaches the disk surface, as illustrated in 
Fig. [6] 

The rather short time delays measured by Sergeev et al. 
(2005) strongly favored disk/corona geometry and low incli- 
nation of the source. High inclination of this source (i = 58^), 
suggested by Crummy et al. (2006) on the basis of the reflec- 
tion model of the soft X-ray excess, is ruled out at the 2-cr level. 
Both models in their simple versions cannot account for a de- 
crease of the variability with wavelength since the irradiation, 
needed to reproduce the observed spectra redder than a bare 
disk spectrum, also increases with the wavelength. However, 
our simple approach to the wavelength-dependent variability 
amplitude did not take into account the increase of the repro- 
cessing region with the wavelength which may result in addi- 
tional suppression of the variability in R and I bands. The pos- 
sibility of a starlight contribution alone did not seem to provide 
the solution. 

The disk/corona set of models favored by the delay mea- 
surement of Sergeev et al. (2005) have black hole masses in 
the range from 1O 7 M to 3 x 10 7 M o . Such values are in agree- 
ment with the best mass measurement based on reverberation 
(M = 1.42 + 0.37 x 1O 7 M , Peterson et al. 2004). The accept- 
able accretion rates cover a broad range (see Table@|i. However, 
if we adopt the Peterson et al. value of the mass, the required 
accretion rate is within the range of 0.55 - 0.75 (in dimension- 
less units). This is consistent with the typical values derived 
for NLS1 objects (e.g. Kuraszkiewicz et al. 2000). The other 
parameters of the model, required by fitting the continuum, are 
also reasonable. The parameter y in all acceptable models was 
1 .8, and the fraction / of the bolometric luminosity undergoing 
reprocessing should be between 0.12 (for an accretion rate of 
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0.55) and 0.07 (for an accretion rate 0.75), for the Peterson et 
al. mass. On the other hand, if indeed a small inclination angle 
is favored for this source, the estimate of the mass from rever- 
beration may be strongly biased (see the discussion by Collin 
et al. 2006). 

The corona model is roughly consistent with the ex- 
pectations of the Inverse Compton heated disk/coronae 
and disk/wind scenarios developed in a number of papers 
(Begelman, McKee & Shields 1983, Ostriker et al. 1991, 
Murray et al. 1994; see Woods et al. 1996 for hydrodynami- 
cal simulations). The inner disk develops a static corona with 
negligible thermal outflow while the outer region, where the 
Inverse Compton temperature is larger than the local virial 
temperature, develops a strong wind. The transition region is 
estimated to be at a fraction of the Inverse Compton radius 
(maximum strength of the corona at 0.2 in the simulations 
of Woods et al.). The disk optical/UV spectra obtained in the 
present paper are similar to those determined from more ex- 
act coronal solutions (e.g. Kurpiewski et al. 1997). Assuming 
an Inverse Compton temperature corresponding to Mrk 335 of 
1.3 X 10 7 K, we obtain the predicted outer radius of the corona 
of4x 10 4 Rschw, of the same order as required by our fits to the 
data. 

The disk/warm absorber geometry is not consistent with 
the delay measurements of Sergeev et al. (2005). This indi- 
cates that there is no massive outflow of the fully ionized ma- 
terial along our line of sight to Mrk 335, although such an 
outflow of considerable optical depth was suggested to take 
place in another NLS1 galaxy (PG121 1 + 143, King & Pounds 
2003). Our results are consistent with the values obtained for 
the partially ionized warm absorber column in Mrk 335 (small 
value of the column density Nh = 4.2+JH x 10 22 crrT 2 , equiv- 
alent to r = 0.028, Bianchi et al. 2001; a similar value of 
N H = 2.5 x 10 22 cm -2 was obtained by Longinotti et al. 2006), 
since such a warm absorber would not backscatter a significant 
amount of radiation. Unfortunately, for PG121 1 + 143, no delay 
measurements exist, so we cannot perform consistency tests for 
that source. 
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